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We recalculate the amplitude for photon splitting in a strong magnetic field 
below the pair production threshold, using the worldline path integral variant 
of the Bern-Kosower formalism. Numerical comparison (using programs that 
we have made available for public access on the Internet) shows that the 
results of the recalculation are identical to the earlier calculations of Adler 
and later of Stoneham, and to the recent recalculation by Baier, Milstein, and 
Shaisultanov. 
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Photon splitting in a strong magnetic field is an interesting process, both from a theoret- 
ical viewpoint because of the relatively sophisticated methods needed to do the calculation, 
and because of its potential astrophysical applications. The first calculation to exactly 
include the corrections arising from nonzero photon frequency uj was given by Adler [1], 
who obtained the amplitude as a triple integral that is strongly convergent below the pair 
production threshold at uj = 2m, and who included a numerical evaluation for the special 
case uj = m. Subsequently, the calculation was repeated by Stoneham [2] using a different 
method, leading to a different expression as a triple integral, that has never been compared 
to the formula of Ref. [1] either analytically or numerically. Recently, a new calculation 
has been published by Mentzel, Berg, and Wunner [3] in the form of a triple infinite sum, 
and numerical evaluation of their formula by Wunner, Sang, and Berg [4] claims photon 
splitting rates roughly four orders of magnitude larger than those found in Ref. [1]. Since 
this result, if correct, would have important astrophysical implications, a recalculation by 
an independent method seems in order. We report the results of such a recalculation here, 
together with a numerical comparison of the resulting amplitude with those of Adler and of 
Stoneham, as well as with a recent recalculation independently carried out by Baier, Mil- 
stein, and Shaisultanov [5]. The comparison shows that these four independent calculations 
give precisely the same amplitude, showing no evidence of the dramatic energy dependent 
effects claimed in Refs. [3] and [4]. 

Our recalculation of the photon splitting amplitude uses a variant of the worldline path 
integral approach to the Bern-Kosower formalism [6, 7, 8, 9]. As is well known, the one 
loop QED effective action induced for the photon field by a spinor loop can be represented 
by the following double path integral, 




(1) 
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Here s is the usual Schwinger proper-time parameter, the x^(t)'s are the periodic functions 
from the circle with circumference s into spacetime, and the ■?/' M (r)'s are antiperiodic and 
Grassmann-valued. 

Photon scattering amplitudes are obtained by specializing the background to a sum of 
plane waves with definite polarizations. Both path integrals are then evaluated by one- 
dimensional perturbation theory, i.e. one obtains an integral representation for the N- 
photon amplitude by Wick-contracting N "photon vertex operators" 

V= f dT[x l "e^-2i^ip ,/ k l ,e v }exp[ikx(T)}. (2) 

J 

The appropriate one-dimensional propagators are 



(^(ri)^(r 2 )) = -g^Gsin, r 2 ) = -gT 



n - t 2 



WfaWfo)) = lg^G F (T U T 2 ) = ±<r*s*(Ti - T 2 ) . (3) 

The bosonic Wick contraction is actually carried out in the relative coordinate y(r) = 
x(t) — x of the closed loop, while the (ordinary) integration over the average position 
x o = 7 Jo drxij) yields energy-momentum conservation. 

To take the additional constant magnetic background field B into account, one chooses 
Fock-Schwinger gauge, where its contribution to the worldline Lagrangian becomes 

AC = \iey^F^ - ie^F^ . (4) 

Being bilinear, those terms can be simply absorbed into the kinetic part of the Lagrangian 

[9 ,10]. This leads to generalized worldline propagators defined by 

1 / d 2 d\ 1 

2ieFi-)g B (T 1 ,T 2 )=5(T 1 -T 2 )-- (5) 



2\dr 2 drj 

\(^-2ieF)G F {T 1 ,T 2 ) = 5{T 1 -T 2 ). (6) 

The solutions to these equations can be written in the form [11] 

Qb{tut 2 ) = -±- 2 ( . y— e~ iesF( ^ B12 4- ieFGBi2 - -\ (7) 
2(eF) Vsm(esF) sj 

„— iesFGnvz 

Gf{ti, t 2 ) = G F12 — — , (8) 

cos(esF) 
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(we have abbreviated Gsij '■= GB(Ti,Tj), and a dot always denotes a derivative with respect 
to the first variable). Those expressions should be understood as power series in the field 
strength matrix. To obtain the photon splitting amplitude, we will use them for the Wick 
contraction of three vertex operators Vq and Vi,2, representing the incoming and the two 
outgoing photons. 

The calculation is greatly simplified by the special kinematics of this process. Energy- 
momentum conservation k + k\ + ki = forces collinearity of all three four-momenta, so 
that, writing — k = k = urn, 

k\ = —k, k<i = —k; k 2 = k\ = k\ = k ■ k\ = k ■ k 2 = k\ ■ k 2 = 0. (9) 

Moreover, a simple CP-invariance argument together with an analysis of dispersive effects 
[1] shows that there is only one allowed polarization case. This is the one where the incoming 
photon is polarized parallel to the plane containing the external field and the direction of 
propagation, and both outgoing ones are polarized perpendicular to this plane. This choice 
of polarizations leads to the further vanishing relations 

£i,2 • £o = £i,2 • k = £i, 2 • F = . (10) 

In particular, we cannot Lorentz contract S\ with anything but e 2 - This leaves us with only 
a small number of nonvanishing Wick contractions: 
1 2 



(V0V1V2) =iexp - LOiUJjnQsijn 
- ij=0 
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^1^512^2 + £l^F12£2^1^2^^F12^ 

x [- ^2 ^i £ oGBOin + iOQe Q F oon - G]&\G] 2 £\Qf\2£2 nQ Fl0 e nQ F20 n - (1 <-> 2) \. (11) 

i=0 J 

For compact notation we have defined u>o = 00,^1,2 = —^1,2- This result has still to be 
multiplied by an overall factor of [^"f^nhlesB) ' wn i c h ^y itself would just produce the Euler- 
Heisenberg Lagrangian, and here appears as the product of the two free Gaussian path 
integrals [8]. 

It is then a matter of simple algebra to obtain the following representation for the matrix 
element C 2 [uj,uji,uj 2 , B] appearing in Eq. (25) of [1]: 
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2eB smh(esB) 
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'cosh(es£>) 

Here translation invariance in r has been used to set the position r of the incom- 
ing photon equal to s. Coincidence limits have to be treated according to the rules 
G b (t,t) = 0,G|(t,t) = 1. 

Alternatively, one may remove Gbi2 by partial integration on the circle. This leads to 
the equivalent formula 

cosh(es-B) 

\2~- 



in r°° 

C 2 [uj , ui, u 2 , B] — — dsse 
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(13) 



cosh(es-B) ^ cosh(es.B) 
This form of the amplitude is less compact, but the integrand (apart from the exponential) 
is homogeneous in the lo^. 

Finally, let us remark that the analogous expression for scalar QED would be obtained 
by deleting all terms in Eq. (11) containing &Qf, as well as the cosh (esB) appearing in the 
overall factor and the global factor of —2 in Eq. (1). 



In order to compare the amplitudes of Eqs. (12) and (13) to those of Refs. [1], [2], and 
[5], we observe that both Eq. (12) and Eq. (13) can be written in the form 

C 2 [u>, ui, u> 2 , B] = —— / — e~ m s J 2 {s, u, ui, u 2 , B) , (14) 

AB z ujujiuj2 Jo s 

in which J 2 is independent of the electron mass m. Inspection shows that the amplitude 
expressions of Adler [1] and Baier, Milstein, and Shaisultanov [5] are already in the form 
of Eq. (14), while that of Stoneham [2] can be put in this form by doing an integration by 
parts in the proper time parameter s, using the identity 

mV m2s = _A e -™ 2 * (15) 
ds 

to eliminate a term proportional to m 2 in the amplitude. In rewriting Stoneham's formulas 
in this form, we note that his Mi(B) is what we are calling C 2 [cu, u>i, cu 2 , B], and that there 
is an error of an overall minus sign in either his Eq. (37) or the first line of his Eq. (40). 
Similarly, in rewriting the formulas of Baier, Milstein, and Shaisultanov in this form, we 
note that their amplitude T is related to C 2 by 

i 8 

C 2 [u>,u 1 ,cj 2 ,B] = T1 ^ T . (16) 

Once all amplitudes are put in the form of Eq. (14), we can compare them by comparing 
the proper time integrand J 2 (s, cu, u>i, u> 2 , B), which in each case involves only a double 
integral over a bounded domain. The only remaining subtlety is that we must remember 
that J 2 vanishes as ujujiuj 2 for small photon energy; this is manifest in Eq. (13) above, but 
in Eq. (12) and the corresponding equations obtained from Refs. [1], [2], and [5], there is 
an apparent linear term in the frequencies which vanishes when the double integral is done 
exactly. In order to get robust results for small photon frequency when the double integral 
is done numerically, this linear term must first be subtracted away, by replacing expressions 
of the form 

"e Q {L + C) (17a) 
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with L, Q, and C respectively linear, quadratic, and cubic in the photon frequencies, by the 
subtracted expression 



This subtraction is already present in the expression of Eq. (25) of Ref. [1], and is discussed 
in the form of Eqs. (17a, b) in Ref. [5], and it also must be applied to Eqs. (37) and (39) of 
Ref. [2] after the integration by parts of Eq. (15) has been carried out. While in principle 
this subtraction should be applied to Eq. (12) above, it turns out not to be needed there, 
because the linear term in the frequencies involves only integrals of the general form 



which is exactly zero using a discrete trapezoidal integration method when the 5 function is 
discretized as a Kronecker delta. Thus Eq. (12) is robust for small photon frequencies as it 
stands, when used in conjunction with trapezoidal integration. 

With these preliminaries out of the way, it is then completely straightforward to pro- 
gram the functions ^(s, ^2> B) for the five cases represented by the formulas of Adler 
[2], Stoneham [3], Eq. (12) of this paper, Eq. (13) of this paper, and Baier, Milstein, and 
Shaisultanov [5], with the result that they are all seen to be precisely the same; the resid- 
ual errors approach zero quadratically as the integration mesh spacing approaches zero, as 
expected for trapezoidal integration. We have not carried out the s and uj\ integrals needed 
to get the photon splitting absorption coefficient, since this was done in Ref. [1], with re- 
sults confirmed by the more extensive numerical analysis given in Ref. [5]. However, anyone 
wishing to do this further computation can obtain our programs for calculating the proper 
time integrand J2 by accessing S. L. A.'s home page on the Institute for Advanced Study 
web site ( [http:/ / www.sns.ias.edu/~adler/Html/photonsplit.htmll ). 




(17b) 




(18) 
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